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We study the melting of a three stranded DNA by using real space renormalization group and 
exact recursion relations. The prediction of an unusual Efimov-analog three chain bound state, that 
appears at the critical melting of a two chain DNA, is corroborated from the zeros of the partition 
function. The distribution of the zeros have been studied in detail for various situations. We show 
that the Efimov DNA can occur even if the three chain [i. e., three monomer) interaction is repulsive 
in nature. In higher dimensions, a striking result that emerged in this repulsive zone is a continuous 
transition from the critical state to the Efimov DNA. 



I. INTRODUCTION 

In recent times the formation of triple helical DNA has 
been a topic of considerable importance because of possi- 
ble implications in the field of molecular biology. In 1957, 
it was discovered that certain sequences of Watson- Crick- 
double helical DNA allow a third strand DNA to bind 
via Hoogsteen or reverse Hoogsteen base pairing to form 
a triple helix[TH3]. This triple helix formation has the 
potentiality to block transcription and thereby affecting 
gene expression. Following this discovery the experimen- 
tal demonstration of the ability of a third chain to recog- 
nize the base sequences without the double helical DNA 
revealing the base pairs, renewed the interest in triple he- 
lix DNA, especially its therapeutic applications [H [5]. It 
is now known that not only DNA, even RNA[B] or PNA 
(polypeptide nucleic acid) are capable of forming triple 
helix with duplex DNA0E]. 

A three stranded DNA has been shown to exhibit an 
Efimov like bound state near or at the critical melting of 
a duplex DNA 9 . The Efimov effect is the most strik- 
ing phenomenon to occur in quantum three-body systems 
with only two-body short-range pair interactions |10IfT5] . 
An infinite number of bound states appear at the critical 
threshold of two body binding. There are several the- 
oretical and experimental investigations using different 
models and methods to show this effect P3HE] . The uni- 
versality of this phenomenon encompasses the analogous 
classical model, namely the melting of three stranded 
DNA[9]. The analogy is drawn between the large quan- 
tum fluctuations near the zero energy threshold of two 
body binding and the thermal Gaussian fluctuations at 
the melting of a duplex DNA. A scaling argument was 
used in Ref[H] to justify the occurrence of the effective 
two chain attractive potential \ as a source of the Efi- 
mov effect. Such a long range interaction leads to a broad 
three-strand DNA bound state at or beyond the melting 
point of a duplex DNA. This is a state where no two are 
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bound but the three are bound together. Such a loosely 
bound state we called the Efimov-DNA. This has also 
been observed from the RG flows and the exact numer- 
ical calculations for several model systems, in particular 
on hierarchical lattices. 

Hierarchical lattices, by virtue of their discrete scal- 
ing, allow one to solve many models in statisti- 
cal mechanics by exact renormalization group (RG) 
transformations [TSHU]- Furthermore many approximate 
real space RG on real lattices can be viewed as exact real 
space RG on hierarchical lattices. In the first study of 
the Efimov effect for Gaussian polymers, RG and exact 
numerics were used [5]. A part of our aim here is to ana- 
lyze the Efimov phenomena exhibited by triple stranded 
DNA from the classical phase transition point of view, 
especially by looking at the zeros of the partition func- 
tions. 

Finding the zeros of a partition function in the com- 
plex plane of any physical variable is a mathematical way 
to understand and analyze the phase transition phenom- 
ena. However finding those is often possible only for 
small sizes or soluble cases and not in general. Yang 
and Lee first studied the Ising Ferromagnetic system in a 
complex magnetic field to show that for a properly cho- 
sen variable the zeros lie on a unit circle, known as the 
Yang-Lee circle [221 125] . Later the zeros were studied in 
the complex temperature plane and other variables |24j. 
Since there cannot be any real zero, the zeros may ac- 
cumulate and pinch the real axis as a limit point in the 
thermodynamic limit. This limit point then identifies a 
transition point. This method can provide relevant in- 
formation on phase transitions such as the critical field 
or temperature and the values of the associated critical 
exponents. Moreover, the distribution of zeros may form 
many complicated structures other than a circle. These 
structures are the separatrices of the two types of flows 
to the two different stable fixed points of the RG trans- 
formation, and are similar to the Julia sets (see appendix 

01231251. 

In Ref[H] the RG flows were studied in the unbound 
region of the two and the three chain state. By looking 
at the flows in the unbound region of the duplex DNA, 
where chains are supposed to be free, an effective three 
chain bound state was predicted. In this paper we study 
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the partition function of the three chain system by com- 
bining the recursion relations and the RG transforma- 
tions, and then finding the zeros. We also extend the 
model to the three chain repulsive interaction regime. In 
addition, we discuss several other features of the zeros in 
the complex plane, for instance the detail structure, and 
the connection to the Julia set. 

This paper is organized as follows. In section [TTJ the 
three polymer problem on a hierarchical lattice has been 
introduced. In section |III[ the recursion relations from 
RG decimation and those for exact iterations are writ- 
ten. The method of finding the zeros of the partition 
function is discussed and we find the limit point of the 
zeros to locate the phase transition. Section |TV| contains 
the results and discussions on the two and the three chain 
system under different situations. In particular we esti- 
mate the transition point for Efimov-DNA. Section |V|ex- 
tends the problem to three chains repulsive interaction. 
The existence of a transition between the Efimov DNA 
and the critical repulsive state in higher dimensions is 
established there. Appendices [A] and [B] are for Julia set 
and limit cycle. 

II. MODEL 

Let us consider the diamond hierarchical lattice as 
shown in Fig. [I] The lattice is generated iteratively by 
the replacement of each bond at the (n — l) th generation 
by a motif of Xb bonds to get the n th generation, where A 
and b represent the bond scaling factor and the branching 
factor respectively. The thermodynamic limit is obtained 
by n — ¥ oo and in that limit the effective dimensionality 
of the lattice is 



In this paper we shall choose A = 2. 



(a) 




n= n= l n= 2 



FIG. 1: (a)The recursive construction of the hierarchical lat- 
tice with b = 2 for n = 0,1,2,... generations. The right 
arrows represent the direction of iteration towards larger lat- 
tices. The left arrows represent the direction of decimation 
used in RG. (b) A motif of 2b bonds, where b = 4. 

One major feature about hierarchical lattices is their 
unusual scale invariance property. They have a discrete 
scaling symmetry. That is why an exact implementation 
of the real space RG technique is possible. The deci- 
mation of the n th generation to arrive at the (n — 1) 
generation is precisely what is needed in an RG trans- 
formation. Once the partition function is known, it is 




FIG. 2: Examples of three-chain configurations on a diamond 
motif for b — 4. (a)Polymers do not share any single bond. 
The number of such configurations is b(b — 1)(6 — 2). (b)Two 
polymers share a bond. The energy here is — 2e and the num- 
ber of such configurations is b(b— 1). (c)Three polymers share 
the same bond. The energy is 2(— 3e — £123). The number of 
such configurations is b. 

possible to calculate the free energy, and the other ther- 
modynamic quantities. One may even write down recur- 
sion relations for them. 

We consider three directed polymers on a diamond hi- 
erarchical lattice. Three chains on the diamond hierar- 
chical lattice are stretched from bottom to top, but they 
can wander at intermediate points. The contact energies 
are defined at the bonds only. The polymers are assigned 
attractive potentials — e and —6123 (e, £123 > 0) if a single 
bond is shared by the two and the three polymers respec- 
tively (see Fig. [2]). At each generation, the length of each 
polymer increases by a factor A = 2 so that the length of 
polymers at the n generation is 

L n = 2". (2) 

For the Efimov effect, only pairwise interaction is 
enough. However in an RG procedure it is imperative 
to define the model with both e and £123, because the 
three chain interaction gets generated on a longer scale. 

III. METHOD 
A. Renormalization Group 

In this subsection we summarize the RG transforma- 
tions and the exact recursion relations for the partition 
functions. The two ways of handling the problem are just 
two different ways to look at it. In the RG case, we start 
from a large lattice and remove short scale fluctuations 
by renormalizing the parameters, effectively reducing the 
size of the lattice. In contrast to this idea of thinning out 
the degrees of freedom, in the second method the lattice 
is built generation by generation so that one may study 
the behaviour of any quantity of interest as a function of 
length of the polymers. This is useful in studying phase 
transitions because finite size scaling can then be used to 
explore the nature of the transition. 

We introduce the Boltzmann factors, 

y = exp(/3e), and w = exp(^e 123 ), (3) 

where j3 — 1/fegT, kg being the Boltzmann constant, 
T the temperature. The RG transformations of the two 
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chain and the three chain Boltzmann factors are given by The initial conditions are taken as 



(b-i) + y 2 



(4) 



C = 1,Z = y,Qo = y 3 w- 



(12) 



, (b-l)(b-2)+3(b-l)y 2 +y b w 2 

W = W ' (5) 

where the primed variables on the left hand side represent 
the renormalizcd values. These recursion relations show 
that the three body term is generated even-though we 
start with £123 = 0, i. e., w = 1. As expected the three 
chain interaction does not affect [i. e., renormalize) the 
two-chain interaction. 

For a given y and w, the flows from successive use of 
Eqs. (|4|-([5]) would give us the phases and the nature of 
the transitions. One needs the fixed points for this anal- 
ysis. For the two chain system, the fixed points of y are 
(i) y* — 00 (stable zero temperature fixed point repre- 
senting a bound duplex state), (ii) y* = 1 (stable infinite 
temperature fixed point representing an unbound state), 
and (iii) y* — (b — 1) (unstable fixed point representing 
the two chain melting or critical point). For pure three 
chain interaction (y = 1) the fixed points of w correspond 
to (i)oo (zero temperature), (ii) 1 (infinite temperature), 
(iii)(6 2 — 1) (unstable, three-chain critical point). The 
two chain melting is critical with a diverging length scale 
with exponent |18j 



log A 



and the specific heat exponent 



a = 2-v. 



(6) 



(7) 



At the two chain critical point y c = b — 1, the fixed point 
of w are found to be 



b 2 ± V4 - 246 + 32b 2 - 126 3 + b 4 

w± = 2ib^ir • (8) 

For b = 4, w± = ^ ii^^p, are complex numbers. In the 
range 3 < b < 8.56 no real roots are found from the three 
chain RG relation (Eq. ([5])). These complex roots lead 
to a limit cycle behaviour, which is intimately related to 
the Efimov effect (see appendix [B]). 



B. Exact recursion relations 

With the trace over all configurations the n th genera- 
tion partition functions for single (C n ), double (Z n ) and 
triple (Q n ) chain systems obey the recursion relations 



C n — bC 2 _ 1 , 



6(6 - + 6^_ l3 

6(6-l)(6-2)C^_ 1 

+36(&-l)C2_i^_H 



(9) 
(10) 

(11) 



The average energy and the specific heat are defined as 

(13) 



w dhrQ" 5E n 

E n = — s , and C n = ——. 

ox ox 



where x is the appropriate variable (y or w as the case 
may be). Though the definitions arc different from the 
actual definitions, proportionality factors are not crucial 
here. 

For given y and w, Eqs. ([9|-( 11 ) give the partition func- 
tions for different L n . The average energy and the specific 
heat can be determined for different L n by writing down 
the recursion relations for derivatives of Eqs. (|9|-(|TT|). 



C. Zeros of partition functions Z n , Qn 

If we take w = 1, i. e., no three body interaction, 
then the partition functions are polynomials in y. In 
general, Z n is a polynomial in y of order L n while Q n 
is a multinomial in y and w. These partition functions 
are then completely described by the zeros, albeit com- 
plex. A phase transition is signaled by a real limit point 
of the zeros. However, the rapid growth of the order of 
the polynomials makes it difficult to implement this pro- 
gramme directly. A different representation is used to 
get the zeros pp]. 

By using the RG transformations of y and u>, the recur- 
sion relations from Eqs. ([9|-(11| can be reduced exactly 
to the forms 



Z n (y) - 6 L "Z„_ 1 (y'), (14) 

(15) 



Q n (y,w) = (b L ") 3 / 2 Q n -x(y',w'), 



with y' , w' given by Eqs. Q, ([5]). These relations can be 
verified by direct substitution and, if necessary, by the 
method of induction. 

Since the zeros determine a polynomial completely, the 
two chain partition functions can be written as 



Z n (y) = b L «- 1 l[(y-q j ) 



(16) 



where qj's are the zeros of partition function Z n (y). 
These zeros appear in complex-conjugate pairs. With 
the substitution of Z n (y), Eq. (14 1 becomes 



b^-'Hiy q 3 ) - fc^-i" 1 n (y' - (17) 
J'=l J'=l 

where q^s are the roots of Z n _\. Then the use of Eq. Q 
gives two roots from each factor on the right hand side, 
so that gj's are the solutions of 



(b-i) + y 2 



(18) 
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zeros obtained from Eq. (18) give a fractal like structure 



q ] =±^bq ] -{b-l). (19) 

It clearly shows that if we know the 2™ _1 zeros of 
Zn-i{y), we wm be able to know the 2™ zeros of Z n (y). 
One may start with the roots of Z\ and generate suc- 
cessively the roots of each generation, by just solving a 
quadratic equation. 

Instead of generating all the roots, a random genera- 
tion is more easily implementable. With an initial value 
yo chosen randomly from the two roots of Z\, the new 
roots are determined by Eq. ( 19 1 . If one of them is chosen 
at random and substituted as cjj, the roots for the next 
generation can be found. Thus after the n th iteration, 
the set obtained is basically the zeros in the complex y- 
plane. These roots are nothing but the zeros of partition 
function found from different sizes of the lattice, which 
in this problem would be equivalent to different lengths 
of polymers. The zeros quickly converge and as n — > oo 
we look for the limit point on the real axis. Apart from 
that, the distribution in the complex y-plane itself is of 
interest. This method has been generalized for the three 
chain system. 



IV. BEHAVIOR OF ZEROS: TWO CHAIN AND 
THREE CHAIN SYSTEMS 

A. Two chain system: b — 4 

For different branching factors, fractal like structures 
are obtained from the zeros of the partition functions of 
the two and the three chain systems. We considered only 
b = 4 as a representative of the range where there are no 
real fixed point along the two chain critical line. For 6 = 4 



(a) 
s 





Re y 



FIG. 3: Plot of zeros of Z n (y) in the complex j/-plane for b = 4 
from (a)exact recursion relation for n = 6, (b)RG relation. 
The closest point to the Re(j/) axis approaches y c = 3, the 
two chain melting point in the limit n — > oo, the unstable 
fixed point of Eq. Q . The selected region shown by a box is 
zoomed in Fig. [2k. 

the structure shown in Fig. [3^, is obtained in the complex 
j/-plane from the exact recursion relation Eq. (14). Ex- 
act solutions are possible only upto n 



6 generation 

because of computational hardware limitations. This is 
insufficient, as thermodynamic limit (n — > oo) is needed 
to observe a phase transition. Finding zeros at random 
from the RG relations (Eqs. Q, (p)])) wins over such dif- 
ficulties and hence large lengths can be reached. The 



shown in Fig. |3p. The accessed zero nearest to the real 
axis approaches the two chain transition point y c = 3 for 
large n. Apart from the limit point, the distribution of 
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FIG. 4: Zeros of Z n (y): The inner rectangular box is zoomed 
successively. A self similar structure becomes apparent. Note 
that the zeros are known with high accuracy. 

the zeros in the complex y-plane is also non-trivial. 

The first feature to note is that the zeros do not seem to 
lie on a smooth diffcrentiable curve. The zoomed picture 
of a small cross-section of the structure for the two chain 
system (from Fig. ^p) has been shown in Fig. [4^,. Fur- 
ther the selected regions have been zoomed successively 
and are shown in Figs. [4]b, [3J:. The self similarity of the 
structure is visible. This is an indication of the fractal 
nature of the distribution. Further analysis required for 
a quantitative description is not done here. 

These fractal like structures obtained above are noth- 
ing but the separatrices of the set of RG flows in the com- 
plex plane to the appropriate stable fixed points. These 
separatrices for iterations of any function in the complex 
plane are known as the Julia set (see appendix |A|). The 
sets are obtained after infinite number of iterations of a 
recursive formula by identifying the points that do not 
flow to the stable fixed points. Our method of finding 
the zeros by using the RG relations is in fact equivalent 
to an inverse iteration method, which is more efficient in 
producing such structures. 

In Fig. [5^i the RG flows are shown in the complex y- 
plane for a two chain system. The dotted line (red curve) 
is showing the flow towards the stable fixed point y = 1, 
i. e., the high temperature region, when we start with a 
value from the inner region of the fractal like structure. 
On the other hand, a point from the outskirts of the line 
of zeros flows to the stable fixed point y = oo, which is 
the bound state with zero temperature. The critical point 
being an unstable fixed point does not actually belong to 
the set but, as discussed, is a limit point — in a sense a 
boundary of the set. 

The second feature to note is the 3-like shape near the 
real axis limit point. It is not arbitrary. The angle at the 
limit point in the complex plane is related to the specific 
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heat exponent by[! 



tan (<pv) 



tan (ira) + — — cosec(7ra), (20) 



where <f> is the angle between the tangent of zeros at the 
limit point, to the real axis of y, and A± are the ampli- 
tudes of specific heat on the low and the high y side of 
the transition. Just like the exponents, A-/A+ is a uni- 
versal number for a universality class of transition. For 
the two chain problem, we know that A^/A + — > oo as 
A + = 0. Therefore the angle <j> is given by 



7T 

2z/ 



(21) 



The zeros obtained by the successive iterations of the 
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FIG. 6: (a)Plot of the zeros of Q n (l,w) in the complex w- 
plane for b = 4. The closest point to the real axis approaches 
w c = 15 for large n. There is self similarity in the distribution 
of zeros, but not shown here. (b)The triangles are the zeros. 
The solid (red)hne given by Eq. ( 22 1 passes through them 
with v of Eq. (27 1 and c = w c . 



where q*j, q/s are the zeros of Q n _i(l,io) and Q n (l,w) 
respectively. The distribution of zeros is the Julia set give 
the fractal like structure shown in Fig. [6^i. By choosing 
the zero near to the limit point w c , the nature of the 
distribution can be determined as shown in Fig . |6p by 
the straight line given by Eq. ( 22 ) with v of Eq. [27] and 

C = W a - 



FIG. 5: (a)Plot of zeros of Z n (y) in the complex y-plane. Two 
types of RG flows are shown. The dotted (red) curve starts 
from a point of the inner region and flows to y = 1. The 
dashed (blue) curve starts from a point of the outer region 
and flows to oo. (b)The triangles are the zeros and approach 
the limi t po int y c = 3 at large n. The solid (red) line, given 
by Eq. ( |22[ ) makes an angle (f> with the real axis with v of 
Eq. [6] and c — y c . 

one close to the real axis are shown in Fig. |5j) by the 
triangles. They approach the real axis in a linear fashion 
with an angle <f>, given by the straight line 



Im z = (Re z — c) tan 



2i/ ' 



(22) 



in generic complex z plane with v from Eq. ([6]). Here c 
represents the limit point of the zeros on the real axis. 
The zeros occur in complex conjugate pairs. Therefore 
if we take the mirror image of the distribution of zeros 
about the real axis in Fig. [5}j, the beak of the 3- like shape 
can be obtained. 



B. Three chain system: b — 4 

We have calculated the zeros of Q n (l,w) for a three 
chain system with a pure three chain interaction. By 
considering y = 1 in Eq. ([5]), we get 



w 



, (b 2 -l)+w 2 



b 2 



(23) 



The zeros come from the equation 



q J =±JWq j -(V-l), 



C. Efimov DNA: b = 4 

The idea is to show the Efimov transition point of DNA 
by finding the limit point of zeros on the real y axis. 
Although we consider w = 1, the effective three chain 
interaction develops by renormalization. As a result the 
zeros found from Eqs. |9])-(11| seem to pinch the Re(y) 
axis at a point where none of the pair is bound. The 
exact solutions are shown in Fig. [7^, for n — 6. On a 
more finer scale the zeros are shown in Fig. [TJd. For such 
small lattices the limit point is not accessible, hence an 
extrapolation scheme may be used. The zeros nearest 
to the Rc(y) axis, obtained in different generations (n = 
2, ...,6) are shown in Fig. [7ji by black dots. A straight 
line nicely fits these zeros, and is shown by a solid (red) 
curve. 

The straight line intersects the real axis at y — 2.321. 
This value is the large n extrapolation and can be taken 
as the estimate of the Efimov transition. We may com- 
pare this extrapolated value with the previous RG based 
estimate of yE — 2.32402. Finding the zeros for two 
chain system is easier than for the three chain system. 
Since the three chain equation holds both the variables 
y, w, finding zeros from the three chain RG relation is 
tantamount to generating the full relation for Q n . This 
is because one needs to keep w at all the intermediate 
values of n and then at the end of the desired value of n, 
w is to be set to 1. One sees the difficulty of the Efimov 
physics eventhough w = 1. It is tempting to simplify 
the recursion relation at the cost of some approximation. 
We set w = w' = 1 to get a renormalized y' that would 
describe the three chain system. Such a relation follows 
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2.16 2.24 2.32 2.4 2.48 
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FIG. 7: Plot of zeros in the complex y-plane for 6 — 4. 
(a)Zeros of Q n (y, 1), when n = 6, (b)a finer scale of (a) near 
the real axis, (c)Combined plot of zeros. The bigger black cir- 
cles are the zeros closest to the real axis (i e., with smallest 
imaginary part) obtained from Q n (y,l) for n = 2,..., 6 and 
the solid straight line is a fit to these. The shaded regime 
shows the distribution of zeros (from Eq. 24 1 on which we su- 



perpose the positive quadrant of (a) shown by the small black 
dots. 



occurs at w — we, from the Efimov DNA to the critical 
state of polymer pairs. Fig|8^ shows the distribution of 
zeros of Q n (y c ,w) in the complex w plane. The set of 
these zeros is a Julia set, separating the flows to the stable 
fixed points. The stable fixed point in the inner region 
of the set is given by Eq. ( 26a ) . The zeros near the real 
axis approach w = we linearly, following Eq. (22) with 
c = we and v of Eq. (28) as shown in Fig|8b. A detailed 



discussion is given in the next section, 
(a) 0.2 
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FIG. 8: (a)Zeros of Q n (y c ,w) in the complex ui-plane for 
6 = 9. The stable fixed point w — w s is sh own by a black 
circle. (b)The solid (red) line is given by Eq. (22 I and passes 
through the zeros shown by the triangles, with v from Eq. ( 28 1 
and c = We- 



from Eq. as 

,3 (6-l)(6-2)+3(6-l)y 2 +y 6 



b 2 



(24) 



The zeros obtained from Eq. ( 24 ) spread out in a "milky 
way" over a region in the complex plane of y. The spread 
makes it difficult to make an estimate of the real axis limit 
point, but one may use the width to put a bound on the 
Efimov transition point (see Fig. [7b) . 



D. Efimov DNA at y c = b - 1: 6 = 9 



V. EFIMOV DNA: RG FLOW AND 
NUMERICAL EVIDENCE 

To explore the robustness of the Efimov effect, we now 
include a three chain repulsive interaction along with the 
pairwise attractive one. The three chain interaction is 
attractive when w > 1 and repulsive for < w < 1. For 
w = 0, representing the hard core three chain repulsive 
interaction, three chains can never be on the same bond 
in this model. 



The study along the critical threshold of the two chain 
melting is quite interesting. No real fixed point for w 
exists for Eq. ^ when b is in the range 3 < b < 8.56 
along the y = y c line. For y = y c , the single parameter 
RG relation is 

, _ (6-2) + 3(6-l) 2 + (6-l) 5 W 2 
W ~ 6^(6-1)2 ' {2b > 

The two fixed points for this case are given by Eq. 
For 6 = 9, these are 

w = w s = 0.0655347.. (stable), (26a) 



w = w E = 0.0926684.. (unstable). 



(26b) 



The unstable fixed point as the phase transition point, 
determines the limit point of the zeros of the partition 
function on the real axis. Hence it can be predicted that 
at the two chain melting point, by tuning w, a transition 



1. 6 = 4 

For 6 = 4 the RG phase diagram is shown in Fig. [9^,. 
The solid (red) line is the separatrix connecting the pure 
three chain transition point (l,w c ) to an Efimov tran- 
sition point for w — 0. Each point on the solid line 
represents the Efimov transition point. In other words 
keeping w fixed, by changing y, we would see a melting of 
a loosely bound Efimov DNA with no pairwise binding. 
The region enclosed between this separatrix (solid red 
line) and the y c = 3 line is the Efimov region and (y, w) 
flows to (l,oo). Below the solid (red) line is the high 
temperature zone of denatured DNA, where RG flows 
are to (1, 1). The region right to the y c = 3 line is the 
two chain bound state. The area below the dashed curve 
where the RG flow takes w to zero when two chain pairs 
are strongly bound, represents a different state where one 
finds a three chain bound state but with no three chain 
contact. The dashed line is then a crossover line, ft re- 
mains to be seen if under some conditions this crossover 
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FIG. 9: RG phase diagram in the y-w plane. (a)For b=4. The 
solid (red) curve and the dashed curve represent the separatri- 
ces, where w flows to two different fixed points on either side 
of the separator. For w = 1, y E = 2.32402 is the Efimov-DNA 
transition point. The filled circles are the Efimov transition 
points for w = 0.5, w = 0.2, w — respectively, obtained 
from Fig. 10 1. (b)For b=9. Along y c = 8, there are two real 
fixed points given by Eq. (|26a[). Solid (red) and dashed lines 



are the separatrices. The filled circles are the Efimov transi- 
tion points for w = 0.15 and w = 0.12 respectively, obtained 
from Fig. |10b. 



line becomes a true phase transition line. 




FIG. 10: The three chain average energy per monomer with 
y from direct computation, (a) For b = 4, the average energy 
curves (marked 1, 2, 3) with the fixed values w = 0.5, w = 0.2, 
w — show first order transitions. (b)For 6 = 9, the average 
energy curves (marked 1, 2) with the fixed values w = 0.15, 
w = 0.12 show first order transitions. Curves (marked 3, 4, 
5) with the constant values w = 0.09, w = 0.065, w = show 
a continuous transition at y c = 8. 



2. b = 9, y c = b - 1 

The RG phase diagram is shown in the Fig. [9Jd for 
6 = 9. In the diagram two separatrices (solid and dashed) 
meet at an unstable fixed point. The two fixed points 
w = w s and w — we are shown in Fig. J9fc>. The presence 
of any unstable fixed point reflects a continuous tran- 
sition along the two chain critical line. Hence we can 
say that by tuning the three chain repulsive interaction 
parameter or temperature in the repulsive zone a tran- 
sition can be induced in the Efimov DNA at the critical 
threshold of duplex binding. The transition is from the 
Efimov state to the critical state of pairs dominated by 
the three chain repulsion. The Efimov region is now re- 
stricted by a separatrix connecting the two unstable fixed 



points (1, w c ) and (y c , we) and the critical line y c = 6—1. 

On the critical line at both the fixed points w = w s 
and w — we, y is a relevant variable (unstable in the y 
direction) . But y does not couple to w in the RG equation 
(Eq. Q). The melting for w < we would be similar to 
the pure two chain melting described by Eqs. (JgJ) , |7]). In 
the y-w plane, (i/ c ,we) is a multicritical point where the 
line of first order transitions goes over to a line of critical 
points. 



3. Data collapse 

We now provide numerical evidences for the above RG- 
based inferences. The exact numerical calculations of the 
average energy and the specific heat by iterating the par- 
tition functions and their higher derivatives for lattices 
of various sizes are done for different fixed values of w. 
lOk for 6 = 4 shows that at 



Fig. 10 1 for 6 = 4 shows that at w = 0.5, w = 0.2 
and w = 0, there are first order transitions. The transi- 
tion points estimated from the point of discontinuity are 
shown by the filled circles in Fig. [9^,. They are on the 
separatrix, and are the Efimov transition points for the 
corresponding values of w. 

The energy curves in Fig. [Top for 6 = 9 with w = 0.15 
and w = 0.12, show first order transitions. These transi- 
tion points are shown by the filled circles in Fig. |9Jd. In 
contrast, the energy curves (marked 3, 4, 5) show con- 
tinuous transitions for w = 0.09, w = 0.065, w = 
respectively at y c — 8. This is consistent with the RG 
prediction of Fig. [9Jd. 

The energy and the specific heat curves are shown in 
Figs. [TTJi, [nja for b = 9, y = 1 and in Figs. [l2ji, [T2J3 
for b = 9, y c = b — 1. Also corresponding finite size 
scaling are shown in Figs. [Tip , [TTJi for b = 9, y = 1 
and in Figs. [12b, [T2"] A for 6 = 9, y c = b — 1. Finite size 
scaling behavior of different thermodynamic quantities 
are described by the length scale exponents. In analogy 
with the Eq. ^ , the exponent to describe the three chain 
transition for y — 1 and y c = 6 — 1 at appropriate critical 
points are given by, 



log 2 



log 



2(b 2 -l) ' 
b 2 

log 2 



(27) 
(28) 



Around a critical point one should see a finite size scal- 
ing. Therefore the average energy and the specific heat 
obeying the finite size scaling can be written in the forms 



E - /^/(i^h-w*!), 
C - L 2/l/ /(i 1/,y |w-u;*|), 



(29) 
(30) 



with appropriate v and w* . In Figs. |llb , |llp l we see 
that the average energy and the specific heat scaled as 
EnLn 1 ^ and C n L n 2 ^ u respectively, when plotted with 
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FIG. 11: (a)For b=4. (a)The three chain average energy 
per monomer versus the corresponding Boltzmann factor for 
chain length upto 2 26 when y = 1. The average energy shows 
a continuous transition at the w = w c . (b)The three chain 
specific heat (C„) per monomer with the corresponding Boltz- 
mann factor. (c)Data collapse of energy. (d)Data collapse of 
specific heat. 
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FIG. 12: For b=9. (a)The three chain average energy per 
monomer versus the corresponding Boltzmann factor for chain 
length upto 2 26 when y c =6 — 1. The average energy shows a 
continuous transition at w = we- (b)The three chain specific 
heat (C n ) per monomer with the corresponding Boltzmann 
factor. Length dependence is shown in the inset. (c)Data 
collapse of energy. (d)Data collapse of specific heat. 



\(w — w*)\Ln U with the v of Eq. (|27j) and w 



w c for 

y = 1, all the data collapse on a single curve for different 
lengths of polymers, where n = 6, 7, , 26. 



Figs. [L2p , |12p show similar plots for the critical line 
{lie = b — 1) with v of Eq. (28) and w* — we- Since the 
specific heat diverges with increasing length, data col- 
lapse is good for the case y = 1. The data collapse for 
the case y c = b — 1 is not so good due to a smoother 
behavior of specific heat at the critical point. These es- 
tablish the weak criticality at w — we- 



VI. SUMMARY 

To summarize, the RG relations and exact recursion 
relations are used to study the three chain system on a 
diamond hierarchical lattice. Our emphasis is on the Efi- 
mov like state exhibited by the three chain system at or 
beyond the two chain melting, where no two are bound, 
and the nature of transitions. Fractal like structures are 
obtained for the zeros of the partition functions. These 
zeros when they pinch the real axis, determine the phase 
transition points. We find that all the transition points 
obtained from RG flows, are in good agreement with the 
zeros of partition function on the real axis. The Efimov 
transition point thus found strengthens the prediction of 
the Efimov like phenomena for the three chain system. 
We have shown that the Efimov effect is exhibited by 
a three chain system even if there is a repulsive three 
chain interaction. A transition can be induced in higher 
dimensions from the Efimov state to the three chain crit- 
ical repulsive state at the melting of a duplex DNA. The 
transition to this three chain critical repulsive state is 
continuous and obeys the finite size scaling law with ex- 
ponents obtained from RG. In the (y,w) phase diagram, 
(Vci w e) is a multicritical point. 

We await experimental evidences for the existence of 
the Efimov DNA. Calorimetric experiments looking for 
the thermodynamic signatures of the Efimov transition 
would be as effective as finding the structure of DNA. 
Again the existence of such a state remains a challenge 
for molecular dynamics and Monte Carlo simulations. 
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Appendix A: Julia set 

The standard definition of a Julia set is the set of points 
on the complex plane which flow to a fixed point (no 
divergence) after a function, e. g., 



c, 



(Al) 



is repeatedly applied, where c is any arbitrary constant, 
and could be real or complex. Let us choose c — 0. The 
fixed point solution for c = are z = 0, 1, oo, where z = 1 
is the unstable fixed point. Here for n — > oo, z n+ \ — > 
0, when starting with \zq\ < 1 and z n +i — > oo, when 
starting with \zq\ > 1. Therefore the unit circle \z\ = 1 
is the boundary between the two stable fixed points z = 
0, oo. The unstable point lies on this boundary. 
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Appendix B: Limit Cycle 

For two successive generations Eqs. ^ will be 

w n - w n+ i = f(w n+ i) - w n+1 . (Bl) 



But if the continuum limit is taken, Eq. (Bl) can be 
written as 



dw 



-{w — w + ){w — W-), 



(B2) 



at the critical line y c = b— 1, where I = InL and L = 2". 
For complex w± — a ± i(3, the solution of Eq. (B2) is 
then 



w = a — f3 tan (3 (In I + < 



(B3) 



where 9 is the integration constant. The above equation 
reflects the periodicity of w in In I with the property 



w(l) = w{l\), where In A 



0" 



(B4) 



Here as I increases w approaches ±oo. This behavior can 
be mapped into a limit cycle in the complex plane with 
a phase factor defined by the equation 



(B5) 



With the help of Eq. (B2| and its derivative, <fi will be 



(B6) 



where </>o is the integration constant. 

Our model on the hierarchical lattice is a discrete 
model. Certainly a limit cycle is obtainable from the RG 
relations in the continuum limit, but it is not straight 
forward to do so in the discrete case. 
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